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Single crystals of (Cd\- x Zn x ) 3 As 2 were synthesized from high-temperature solutions and char¬ 
acterized in terms of their structural and electrical properties. Based on the measurements of 
resistivity and Hall signals, we revealed a chemical-doping-controlled transition from a three dimen¬ 
sional Dirac semimetal to a semiconductor with a critical point x„ ~ 0.38. We observed structural 
transitions from a body-center tetragonal phase to a primary tetragonal phase then back to a body- 
center tetragonal phase in the solid solutions as well, which are irrelevant to the topological phase 
transition. This continuously tunable system controlled by chemical doping provides a platform for 
investigating the topological quantum phase transition of 3D Dirac electrons. 

PACS numbers: 06.60.-c, 61.05.-a, 71.30.+h, 72.90.+y 


I. INTRODUCTION 

Cadmium and zinc pnictides ( Cd 3 As2 , Zn 3 As2, Cd 3 P-2 and Zn 3 P2) belong to a group of A 3 r B 2 semiconductors and 
semimetals, which have been well known for their potential applications in high efficiency solar cells and optoelectronic 
devices 14 . These four compounds crystallize at various temperatures in several closely related structures, which can 
be viewed as the different arrangements of a distorted antifluorite structural 0 } 

The electrical properties of these four compounds are distinct in several aspects. Zn 3 As2, Zn 3 P2 and Cd 3 P2 are 
semiconductors with low carrier mobility and direct band gaps being 1.0 e V, 1.5 eV and 0.5 eV respectively 1113 . 
Both Zn 3 P 2 and Zn 3 As 2 are p-type, while Cd 3 P 2 and Cd 3 As 2 are n-type 11 '"™. On the other hand, previous studies 
on the optical properties of Cd 3 As2 suggested that it was a semiconductor with a narrow band gap around 0.1eV^4 
The mobility for Cd 3 As2 was reported as high as 1.5 x 10 4 cm?/Vs at room temperature 11 . For comparing, the hole 
mobility for Zn 3 As2 was only 10 cm 2 /Vs at room temperature. Cd 3 As2 was believed to manifest an inverted band 
structure due to the spin-orbital coupling (SOC) while the other three had normal band structured^. 

Very recent studies on Cd 3 As2 have revealed the topological aspect of its electrical propertied®®!. B an d structure 
calculation predicted that Cd 3 As2 was a three dimensional (3D) Dirac semimetal with the band inversion 15 . The 
energy dispersion of the Dirac electron is protected by the rotational symmetry along the crystallographic c axis 
in the tetragonal unit cell. The 3D Dirac cones of Cd 3 As2 have been observed in angle-resolved photoemission 
spectroscopy (ARPESj^SU. The Dirac-like band dispersion and the inverted band ordering was probed by the 
Landau level spectroscopy and quasiparticle interference in scanning tunneling microscopy (STM) 19 . Based on the 
electrical transport measurements, two experimental groups found an ultrahigh mobility of the Dirac electrons- 0 2 4 
A strongly sample-dependent, large linear magnetoresistance (MR) was observed in Cd 3 As2 at low temperature^!. 
The nonsaturating linear MR in n-type Cd 3 As2 up to 65T was believed to result from its mobility fluctuations^. 
The anisotropic Fermi surface with two elli pso ids of Dirac electrons along the c axis was revealed from the angular 
dependent measurements of SdH oscillations-4 

Noticing the opposite band orderings in Cd 3 As2 and the other three members in the family, we expected a band 
inversion transition in pseudo-binary compounds of (Cd\- x Zn x ) 3 As2 and Cd 3 (As\- x P x )2- A band inversion transi¬ 
tion due to the change of the SOC str ength has been observed in solid solutions of semiconductors and semimetals 
such as Hg 1 _ x Cd x Te and Pbi_ x Sn x S d 2 A 25 [ Recent studies on the solid solutions of TlBiSe2~ x S x and Bi2- x In x Se 3 
have confirmed the e xisten ce of a quantum phase transition tuned by chemical doping from topological insulators to 
trivial band insulatoriP^ 27 . A topological phase transition from a 3D Dirac semimetal to a trivial semiconductor was 
predicted in Na 3 Bii_ x Sb x and Cd 3 (As\- x P x )2 by first-principle calculations as well 28 . 

The changes of the structures and physical properties of polycrystalline (Cdi_ x Zn x ) 3 As2 and Cd 3 (Asx- x P x )i 
have been studiecP^H. (Cd\- x Zn x ) 3 As2 crystallize in a primar y tet ragonal structural. The majorities undergo a 
crossover from n- to p-type when x increases in (Cd\- x Zn x ) 3 As 'P G ^ il \ while the band gap increases linearly with the 
proportion of Zn, according to the magneto-optical measurements^. Due to the lack of single crystals, the change of 
the band structure and the electrons, especially in the topological aspect, cannot be addressed in detail. 

In this study, we report the single-crystalline (Cd\- x Zn x ) 3 As2 obtained from high-temperature solution growth. 
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Powder X-ray Diffraction (XRD) measurements revealed structural transitions from a body-center tetragonal phase 
for Cd 3 As 2 , to a primary tetragonal phase for 0.07 < x < 0.52, and then back to a body-center tetragonal phase for 
x > 0.52. The electrical resistivity and Hall measurements revealed a metal-insulator transition at a critical point 
x c ~ 0.38. The analysis of the MR demonstrated a transition from a 3D Dirac semimetal to a trivial direct-gap 
semiconductor via the modulation of the SOC strength. 


II. EXPERIMENT 

Single crystalline ( Cdi- x Zn x ) 3 As 2 samples were grown from high temperature solutions with the initial concentra¬ 
tion of starting elements being ( Cdi- y Zn y )gAsi. The mixtures were sealed in evacuated quartz ampoules, and then 
kept at a high temperature between 800°C and llOO 0 ^ for two days, and then slowly cooled down to 425 0 C with 
a rate of —5°C/hour. After staying at 425°C for one day, the ampoules were centrifuged to separate crystals from 
flux. The single crystals of Cd 3 As 2 were mainly 3D bulks with triangular facets, but several needle-like crystals were 
found in the growth as well. The shapes of the crystals were similar as that described in Ref. [20,21]. When zinc was 
added to the solutions, the sizes of the 3D crystals decreased, while some flake-like crystals with smooth or mesa- 
landscape-like surfaces appeared. For x > 0.58, the flake-like crystals were dominant and no 3D crystals appeared 
in the growth. XRD measurements revealed that both the triangular facets (Fig. 1(f)) and the large surface of the 
flake-like crystals (Fig. 1(e)) were either the (011)p face of the P4 2 /nmc structure or the (112)/ face of the Mi/cd 
structure, which were the counterparts of each other (Fig. 2(c)-2(d)). In a same batch of growth, the crystals with 
different morphologies did not show larger difference of physical properties than those with the same morphologies. 

In order to determine the zinc concentration x , we measured the Energy Dispersive X-ray Spectrum (EDX) of the 
samples in an FEI Nova NanoSEM 430 spectrometer. The samples with no residual cadmium were selected in the 
measurements and their EDX spectrum was observed through an overall area scanning. The linear relation of the 
measured zinc concentrations x and the initial y (0 < y < 0.1) is shown in Fig. 3(c). By weighing the mass of the 
crystals yielded in every growth, we found that all the initial stoichiometric zinc was compounded in the crystals. 

The powder XRD data was collected from a Rigaku MiniFlex 600 diffractometer and then refined by a Rietica 
Rietveld program. As shown in Fig. 3(b), the lattice constants and the volume of the unit cell change linearly with 
x in accordance with a Vegards law. This result is same with what observed in polycrystalline (Cdi- x Zn x ) 3 A S5 P. 
Therefore we selected x determined by EDX measurements as the nominal zinc concentrations, which were estimated 
to have less than ±1% difference between the real zinc concentrations. More details of the XRD experiments are 
discussed in the Result part. 

Single crystals were polished to bars with length ~ 1.0mm, width ~ 0.4mm and thickness ~ 0.3 mm for elec¬ 
trical transport measurements. The crystallographic orientation of the bars was same as those chosen in the recent 
experimenters), i.e., the current was perpendicular to the [112]/ direction in I\\/cd (or [011]p direction in P4 2 /?/toc), 
and the magnetic field was along the [112]/ direction in IAi/cd (or [011]p direction in P4 2 /nmc). The electrical re¬ 
sistance and Hall voltage were measured via a four-point contact method in Quantum Design Physical Property 
Measurement System (PPMS-9). 


III. RESULT 

The low temperature phases of Cd 3 As 2 and Zn 3 As 2 were reported as a" (P4 2 /n mc), a' (P4 2 /n6c) and a (74/cd) 
at different temperature^' 7 ( which evolve from a high-temperature /3 (Fm3m) phas^H. The /3 phase belongs to an 
antifluorite structure in which an arsenic atom is coordinated by six cationic and two vacancies randomly distributed 
in corners of a cube. In the low-temperature phases, two cationic atoms are missed along a diagonal of one face in 
the distorted cub^A 7 -. At room temperature, both Cd 3 As 2 and Zn 3 As 2 were reported to crystallize in a body-center 
tetragonal phase 6 'EaHSl. Recently single-crystal XRD measurements suggested that the crystals of Cd 3 As 2 form in 
the structure of IAi/acd instead of 74 /cd at room temperature 9 . The unit cell of 74 /cd phase is made of the unit cells 
of P4 2 /nmc phase associated with the lattice constants a/ = y/2ap and c/ = 2cp (The subscripts 7 and P present 
the body-center and primary space group respectively). 

Our powder XRD measurements revealed that {Cd\- x Zn x ) 3 As 2 have different crystal structures for different x at 
room temperature. Figure 2(a) shows that the crystals of Cd 3 As 2 grown from flux are a phase, and their XRD pattern 
exhibits the characteristic peaks of (231)/, (233)/ and (237)/_of the body-center tetragonal phase for 25° <29 < 36°. 
This result is consistent with what was previously reported 7 ®!. Limited by the quality of the data, we could not 
determine whether Cd 3 As 2 crystallizes in the space group of 74/cd or IAi/acd. Once a small amount of Zn was added 
(x = 0.07), the XRD pattern is distinct from that of Cd 3 As 2 ■ The characteristic peaks of (231)/, (233)/ and (237)/ 
disappear, while the (032)p peak of the P4 2 /ranc group occurs (Fig. 2(a)). This peak remains resolvable until the 
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doping level reaches x = 0.38. For 0.38 < x < 0.46, the structure reenters the body center tetr agon al structure IA\cd 
accompanied by (240)/ and (244)/ peaks which are exceedingly weak in the pattern of C d 3 As-p*^(F\g. 2(b)). 

The peaks of (440)/ with the strongest intensity stand at 40.0° and 43.1° for a — Cd 3 As 2 and Zn 3 As 2 respectively 
(Fig. 3(a)). The peak of (040)p is the counterpart of (440)/ in PI 2 /nmc group. Figure 3(a) shows that the peaks of 
(040)p and (440)/ shift gradually when x changes from 0 to 1. Although the volume of the unit cell changes about 
20% from x = 0 to 1, the peak shape does not change significantly in the solid solutions, indicating homogeneous 
chemical distributions in the crystals. The lattice constants for the samples in the 74/ cd group were presented in the 
view of the P4 2 /nmc group as a/ = y/2 ap and c/ = 2cp. Figure 3(b) shows that ap and cp change in a precisely 
linear relation with respect to cc, albeit the structural transitions. 

Temperature dependent resistivity of ( Cdi- x Zn x ) 3 As 2 shows a clear change from a metallic to semiconducting 
profile when x increases from 0 ~ 0.31 to 0.38 ~ 0.58 (Fig. 4(a)). The resistivity of Cd 3 As 2 is close to what was 
reported in Ref. [22] with the residual resistivity ratio ( RRR = P(t=300K)/ P{T=2K)) being 10 (Fig. 8). The values of 
RRR keep almost invariant when x increases up to 0.31 (Fig. 8). For x = 0.38, the p(T) decreases with decreasing 
temperatures above 200AT, and then increases below this temperature. This complicated behavior indicates that the 
sample is likely a very narrow bandgap semiconductor. The values of RRR then dramatically decrease for x > 0.38, 
being 0.58 for x = 0.38 and 7.7 x 10~ 6 for x = 0.58 (see more details in Fig. 8). The changes of the RRR for different 
x indicate a process of band gap opening for x > 0.38. 

Figure 4(b) shows the MR for t he sa mples for x < 0.46 at 2 K. When x < 0.31, the values of the AIR are 
comparably large as that for Cd 3 As 2 R2I251 The va i ues 0 f the AIR decline significantly for the semiconducting samples 
for x > 0.38. Recent studies of Cd 3 As 2 reported linear MR at low temperature d 22 ! 2 ^ . In our experiments, the MR 
of {Cdi- x Zn x ) 3 As 2 follows the power law of MR oc H a where a varies from 0.9 to 1.5 for different samples. 

Conspicuous SdH oscillations occur in the field dependent resistivity for all the samples for x < 0.38 at low 
temperatures, while the oscillations were not observed for the samples for x > 0.46. This observation is consistent 
with a metal-insulator transition with a critical point x c ~ 0.38. The part of resistivity with oscillations versus the 
reciprocal of the magnetic field is presented in Fig. 5. For x < 0.29, only one frequency was observed for each 
sample (Fig. 6(a)). The frequencies of the oscillations show a clear trend of a decline with respect to x up to 0.29 
(Fig. 6(a)). For 0.29 < x < 0.38, the frequencies show more significant sample difference in a same batch. Some 
samples show single frequencies from 15T to 30T, while the second and third frequencies as large as 70 T occur in 
other samples. Such strong sample-dependence and complicated multi-frequency features indicate that the samples 
for 0.29 < x < 0.38 are semimetals or very narrow bandgap semiconductors with complicated Fermi surface which is 
strongly influenced by subtle changes of chemical potential. 

As shown in Fig. 5, the temperature dependent amplitudes of the SdH oscillations were fitted by the Lifshitz- 
Kosevich formulaP^E 


A p xx oc A(T)e 


cos27t(^ + /3) 

JD 


A(T) = 2 t T 2 k B T/hv c 

sinh(2TT 2 kBT/hco c ) 

where kp is the Boltzmann’s constant; w c is the cyclotron frequency; Tp> is the Dingle temperature and A(T) is 
the thermal damping factor which helps to fit the energy gap Hu> c . For the samples with multi-frequencies, their 
main frequencies were analyzed. For large x , the amplitudes of the oscillations are damped less significantly by the 
temperatures for large x , which indicates a smaller cyclotron effective mass (Fig. 5). 

The parameters of the SdH oscillations for different x are listed in Table I. The cross-sectional area Ap in the 
momentum space comes from the Onsager relation Sf = ^;Ap. Simply assuming a circular Ap, we got the Fermi 
wave vector kp from Ap = pk 2 F . The Fermi velocity vf = hkp/m *, the Fermi energy Ep = u F m*, and the cyclotron 
effective mass in* = eB/ uj c are listed as well. Figure 6 shows that m* changes in a similar manner as Sf\ with respect 
to x. 

In order to better understand the metal-insulator transition, we measured Hall resistivity in (Cdi- x Zn x ) 3 As 2 
at 2K. The field dependent Hall resistivity of the samples for x < 0.31 shows a linear negative profile with SdH 
oscillations on the background. The negative linear-field-dependent p yx (H) indicates that the carriers in the samples 
for x < 0.31 simply originate from an electron band. The carrier density decreases linearly with increasing x from 0 to 
0.31 (Fig. 7(a)). These results are consistent with the observation of the decreasing SdH oscillation frequencies with 
respect to x. p yx {H) becomes smaller and nonlinear for 0.38 < x < 0.59. This nonlinear feature is clear for x = 0.46 
(inset of Fig. 7(a)). In this range, the samples manifest a semiconducting p(T ) profile while multi-frequencies were 
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observed in the SdH oscillations in their MR. The Hall signals and the resistivity indicate two types of carriers. 
For x > 0.59, the Hall signals become large and positively Held-dependent, which indicate p-type semiconductors 
consistent with low carrier concentrations. The change of the carrier density rip and mobility p c with respect to x is 
summarized in Fig. 7(b). 


IV. DISCUSSION 

Our measurements show that the change of the electrical properties of {Cd\- x Zn x )zAs 2 has no observable correla¬ 
tion with the structural transitions. This result is not unexpected according to previous band structural calculation. 
Both a” (P^/nmc) and a (J4icd) phases of CdsAs2 manifest similar simple band structures near the Fermi surface. 
Only two Dirac cones protected by rotational symmetry cross the Fermi level (Ep) along the high symmetric line 
r — Z in the Brillouin zonesP^. Therefore a structural transition cannot influence the bands near the Ep. 

Prev ious s tudies showed that the calculation and experiments, the negative gap is about — 0.3eV ~ — 0.7eU for 
Cd^As^^- 1, while the direct gap is l.OeV for Zn^As^-. With a semimetal and a semiconductor as two terminals, the 
band inversion transition should accompany a metal-semiconductor transition at a certain x. If we assume that the 
band gap of {Cdi- x Zn x )^As 2 changes linearly with respect to the critical point of the band inversion transition 
is estimated to occur in the range of 0.23 < x < 0.41. This estimation is consistent with our experimental results. 
The critical point can also be estimated by considering the change of the SOC strength in ( Cd\- x Zn x )^As 2 ■ Here 
we assume that the band inversion is solely induced by the change of the SOC strength, which is proportional to 
Z 4 /n 3 in case of the hydrogenic wavefunctions in a Coulomb field where Z is the nuclear charge and n is the principal 
quantum numbed 39 . Then the critical point is estimated as x ~ 0.35, which is very close to the experimental result: 
x c ~ 0.38. 

The band structure calculation for Cd^,{Asi- x P x )2 revealed a topological phase transition from a Dirac semimetal 
to a trivial semiconductor induced by the change of the SOC strength—1 When x increases, the two Dirac points along 
the k z axis gradually move closer and then merge at the T point under the protection of the crystal symmetry. A 
direct band gap is opened beyond the critical point. The process of the band inversion transition in ( Cdi- x Zn x ) 2 ,As 2 
should be similar as that in C dz{As\- x P x )2- Despite of large chemical replacement in the crystals, the Dirac cones are 
robust, which is supported by the vanishing disorder self-energy around the crossing points 49 . Further investigation 
such as ARPES measurements for (C'd\- x Zn x )^As2 will help to reveal the details of this topological phase transition. 

Cd^A s2 is a lways n-type due to As vacancies, while Zn 3 As2 is p-type because extra Zn vacancies serve as electron 
acceptors®^. Since both two types of carriers come from element vacancies, an n to p transition is expected in 
( Cd\- x Zn x )^As2 ■ With increasing x, the zinc doping will suppress the chemical potential, which crosses a small 
Fermi surface near the Dirac cones. The decrease of Sf i with increasing x is a comprehensive result of the change of 
the band structure and chemical potential. 

For 0.29 < x < 0.38, we found strongly sample-dependent frequencies of SdH oscillations. For a very narrow 
bandgap semiconductor or semimetal, any small change of the carrier concentrations will affect the chemical potential 
dramatically near the band touching. The strong sample-dependence and the complicated SdH oscillations in this 
regime are not unexpected. 


V. SUMMURY 

Single crystals of (Cdi- x Zn x ) 2 ,As 2 were synthesized from high-temperature solutions. Based on the analysis of the 
electrical properties, we realized a transition from a 3D Dirac semimetal to a semiconductor with the critical point 
x c ~ 0.38 in these solid solutions. The structural transitions do not affect the electrical properties in this system. 
The topological aspect of this metal-insulator transition needs experimental exploration in the future. 
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FIG. 1: Panel (a)-(d): Unit cells of a" — Cd, 2 ,As 2 and a — CdzAs 2 viewed along c and b axes respectively. Panel (e) and (f): 
Single crystals of ZnzAs 2 and (Zno. 3 iCdo.e 9 ) 3 As 2 show different morphologies (see more details in the text). The scale on 
right is 1 mm. 
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FIG. 2: XRD patterns for {Cd 1 - x Zn a .) s As 2 - Panel (a): Powder XRD patterns for x = 0, 0.07 and 0.14 for 25° <29 < 36°. 
The characteristic peaks of (231)/, (233)/ and (237)/ for a body-center structure and (032)p for a primary structure are labeled. 
Panel (b): Powder XRD patterns for x = 0.38, 0.52 and 0.58 for 31° <29 < 37°. The peak of (032)p occurs for x = 0.38, while 
the peaks of (240)/ and (244)/ occur for x = 0.58. For x = 0.52, we observed two sets of peaks, which indicates the batch has 
the mixture of two types of crystals. The (100) peak of cadmium at 34.5° occurs in all the XRD patterns with no shift. Panel 
(c): The diffraction pattern of the (112)/ plane for a — Cd^,As 2 and a — Zn^Asi- Panel (d): The diffraction pattern of the 
(022)p plane for {Cdo.sjZno.\z)3As2- 



























Nominal Zn concn. (x) 



Initial Zn concn. (y) 


FIG. 3: Panel (a): The (440)/ peaks and the counterpart (040)p peaks in the powder XRD patterns for (Cd\- x Zn x ) 3 As 2 - 
Panel (b): The lattice constants and the volume of the primitive cell change with the concentration (concn.) of zinc x linearly. 
Panel (c): The nominal concentration of zinc x and the initial y have a linear dependence for 0.00 < x < 0.58. 
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FIG. 4: Panel (a): Temperature dependent resistivity of {Cdi- x Zn x ) 3 ,As 2 for 0 < x < 0.58 in zero magnetic field. Inset: The 
values of resistivity for 0 < x < 0.46 at 2 K in zero magnetic field. Panel (b): MR versus magnetic fields at 2 K. Inset: The 
values of MR for (Cd\- x Zn x )3 As 2 (0 < x < 0.58) versus x at 9 T at 2 K. 
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FIG. 5: The oscillatory components of A p xx versus the reciprocal of the magnetic field (1/B) at different temperatures for 
x = 0.07,0.29 and 0.38. Insets: Fast Fourier Transform (FFT) spectra for x = 0.07,0.29 and 0.38. The results of the SdH 
oscillations are summarized in Table I. 
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Nominal Zn concn. (x) 

FIG. 6: Panel (a): The main frequency of the SdH oscillations at 2I\ for {Cd\- x Zn x )zAs 2 - Different colors represent different 
samples from the same batches. The black ones are the samples corresponding to Fig. 4(a) and 7(b). Panel (b): The cyclotron 
effective mass m* according to the main frequency for each sample versus x. The line is for visual guidance. 
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Nominal Zn concn. (x) 

FIG. 7: Panel (a): Field dependent Hall resistivity of {Cdi- x Zn x )sAs 2 for 0 < x < 0.46 at 2K. From x = 0.00 to 0.46, all 
samples are n-type. Inset: The Hall resistivity of the sample for x = 0.46 shows a non-linear profile. For x = 0.58, the signal 
turns to p-type. Panel (b): Carrier density uh (uh = B/(ep yx )) and mobility p c (p c = l/(ep xx n_fr)) at 2 K versus x for x = 0.00 
to 0.38. The solid line is for visual guidance. 
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FIG. 8: Phase diagram for (Cdi-^Zn^aAs?. With the concentration of Zn increasing, the samples transform from 
a topological Dirac semimetal to a semiconductor. The upper sketches of the band structure illustrate this transition. The 
structure transforms from J4icd to PAi/nmc then back to J4icd with vertical dashed lines serving as rough boundaries. The 
background color presents the gradual change of resistivity at 2 K as x increases (inset of Fig. 4(a)). The changes of Sfi and 
RRR are plotted in the diagram as well. 
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TABLE I: Parameters of the tested samples of different zinc concentration x. The dashed entries mean quantities 
missing. 


'- -X 

Parameters 

0.00 

0.07 

0.13 

0.29 

0.31 

0.38 

A f (*10‘ 3 A' 2 ) 

8.08 

6.18 

2.53 

1.30 | 1.74 

1.98 

2.66 | 4.01 | 6.66 

ma 1 ) 

0.051 

0.044 

0.028 

0.02 | 0.024 

0.025 

0.029 | 0.036 | 0.046 

Ep(eV) 

0.346 

0.289 

0.165 

0.107 | - 

0.143 

0.169 | - | - 

v F (*10 5 m/s) 

9.9 

10.1 

9.1 

8.3 | - 

8.6 

9.4 | - | - 

m* 

0.061 

0.051 

0.036 

0.028 

0.034 

0.036 




